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ABSTRACT 

Context. We present a catalog of 857 white dwarf (WD)-M binaries from the sixth data release (DR6) of the Sloan Digital Sky Survey 
(SDSS), most of which were previously identified. For 636 of them, we complete a spectral analysis and derive the basic parameters 
of their stellar constituents and their distances from Earth. 

Aims. We attempt to measure fundamental parameters of these systems by completing spectral analyses. We propose to test models 
typically applied in fitting procedures and constrain likely and appropriate evolutionary scenarios for the systems. 
Methods. We use a minimization technique to decompose each combined spectrum and derive independent parameter estimates 
for its components. The possibility of alignment by chance is demoted to statistical insignificance, hence, we use physical interaction 
of the binary constituents as input parameter. Additionally, we check the corresponding photometric data from the SDSS to find 
optically resolved systems. 

Results. Forty-one of the stellar duets in our spectroscopic sample are optically resolved in their respective SDSS images. For these 
systems, we also derive a minimum true spatial separation and a lower limit to their orbital periods, typically which are some 10^ yr. 
Spectra of 167 stellar duets show significant hydrogen emission and in most cases no additional He i or He ii features. We also find that 
20 of the 636 WDs are fitted to be DOs, with 16 measured to have around 40 000 K. Furthermore, we identify 70 very low-mass 
objects (VLMOs), which are secondaries of masses smaller than about 0.1 M©, to be candidate substellar companions. 
Conclusions. Although various selection eff'ects may play a role, the fraction 6.4 % of WD-M star binaries with orbital separations 
of around 500 AU is a criterion for evolutionary models of stellar binary systems. Of the 167 spectra with hydrogen emission, 8 had 
already been found to be post-common envelope binaries (PCEBs) and 4 are systems with strong irradiation processes on the M 
dwarf. The remaining 155 B aimer-emitting binaries probably harbor an active M dwarf (dM), corresponding to a fraction of 24.4 %. 
The excess of cool DOs is most likely due to additional WDs in the DB-DO T^^ range, for which no detailed fitting was completed. 
The trend of the M stars being closer to Earth than the WD component is probably due to an underestimation of the theoretical M star 
radii. 

Key words. (Stars:) binaries: spectroscopic - (Stars:) white dwarfs - Stars: late-type - Stars: fundamental parameters - Methods: data 
analysis - (Stars:) binaries: visual 



1. Introduction 



The study of white dwarf (WD)-M star systems is an impor- 
tant research area because they constitute a common final stage 
object of stellar evolution, a WD, and the most frequent type 
of star: the M component. Close WD-M binaries are also pro- 
genitor candidates for cataclysmic variables (CVs) and Type la 
supernovae. WD-M binaries consist of two stars of radically dif- 
ferent structure and evolutionary stage that originate in interstel- 
lar matter (ISM) of identical composition at the same time, on 
cosmological scales, and in the same region of space. Many of 
the constituents of the systems do not only interact gravitation- 
ally, which may lead to mutual mass exchange and substantial 
gravitational radiation, but they also aflfect each other due to 
their magnetic fields. With this paper, we attempt to improve the 
knowledge of the basic parameters of WD-M systems. 



The Sloan Digital Sky Survey (SDSS) provides an invalu- 
able observational data set with which to complete such a study. 
iRavmond elaP (l2003h first attempted to study WD-M dwarf 
(dM) pairs in the SDSS and identified 109 of these objects. 
In our master sample, there are 99 of them but we discarded 
32 for our purposes because of reasons given at the end of 
Sect. [21 ^ilvestri et al. (2 007,) presented a study of 1253 de- 
tached close binary systenis from the SDSS data release 5 (DR5) 
(lAdelman-McCarthy et al.ll2007l) . most of them which consisted 
of WD primaries and M dwarf companions. Our master sam- 
ple comprises 857 WD-M sta r system s, 690 of wh ich have been 
publis hed by ISilvestri et al.l (l2007l) . iReba ssa- Mansergas et al.l 
(l2007l) presented investigations of 101 white dwarf main- 
sequence binaries (WDMSs) with multiple spectra from the 
SDSS, using model atmosphere fitting for the WD primaries 
and spectral typi ng based o n the M-dwarf te mplate spectra of 
iBeuermann et al.l ([l998.) and lBochanski et all 02007.) for the sec- 
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Fig. 1. Left. The localizations of our master sample in the {u - r)-space; Middle: The localizations of our master sample in the 
(/ - z) -space; Right: Comparison of the color-color localizations of our master sample with the WD-M sample of ISilvestri et alJ 
(l20Q6h . Those objects from our study with g - r > 0.2 come from the EDR, on which we did not impose the {g - r)-cutoff. 



ondaries. Significant radial velocity variations were found for 
18 systems, indicative of sufficiently close orbits to be post- 
common envelope binaries (PCEBs), and Sect. 15.21 is devoted 
to this su bject. Another study of WD-M binaries was com- 
pleted by 'Wach ter et al.l (l2003l) . who focused on the observa- 
tional identification methods with data from the Two Micron 
All-Sky Survey (2MASS) increm ental data release . Catal ogs of 
W Ds in the SDSS were gi ven bv iKleinman et al.l (|2004 DRl) 
and lEisenstein et all (l2QQ6l DR4) and many are included in our 
sample. Four binaries from our master sample are in Luyten's 
WD catalogs for the years from 1970 to 1977 (.Luvt en 1999), 
which contains 6546 objects, an d no binary is in his sample of 
6210 objects from 1940 to 1987 (lLuvtenlll997h . The original ac- 
curacy of the position in J1950 coordinates was 0.1 min in Right 
Ascension (1.5' or 0.025°) and \ ' in Declination (0.017°), as 
stated in the respective VizieR reference. We used the J2000 co- 
ordinates calculated by VizieR to relate stars with each other. 
Another earlier stu dy of WDs in wide binaries was completed by 
iGreensteinI (Il986h . and we found three of his 56 systems to be in 
our master sample, one of a quality too low for a reliable analy- 
sis. We also show fits for the three binary systems mentioned in 
iHiigelmeyer et al.l (l2006l) . From our sample, 153 systems were 
not included in any of the mentioned studies. 

The present work introduces a novel method of analyzing 
WD-dM binaries by fitting simultaneously model spectra for 
both the white dwarf and the main- sequence star to the com- 
posite spectrum to derive the atmospheric parameters of both 
stars. Unlike the aforementioned studies, we thus obviate the 
need to calibrate empirical relations between physical parame- 
ters and spectral types. Although current cool-star atmosphere 
models may still have some shortcomings in fitting M-dwarf 
spectra, this approach avoids possible biases due to systematic 
differences between single M-dwarfs and those in close binaries. 

In Sect. [21 we provide an overview of the observational prin- 
ciples and our target selection criteria. We also state the premises 



for our procedure, e.g. we calculate the probability of WD-M 
alignment by chance (Sect. 12. ll) . The calculation of the theoret- 
ical spectra used in the fitting process is described in Sect. [3l In 
Sect.m we then provide a detailed description of the mathemati- 
cal methods we used: This represents the major part of this paper, 
since the method that we use differs from those typically ap- 
plied to decompose combined spectra. With our relatively large 
sample of 636 systems, we present a statistical refinement of the 
SDSS WD and the M dwarf temperature functions in Sect.[5l We 
also find 41 optically resolved binaries with accordingly wide- 
separated constituents of long orbital periods (Sect. 15.31) . 



2. Observations 

The observational data were taken from the SDSS DR6, which 
was performed using the 2. 5 -meter telescope at the Apache 
Point Observatory in southern N ew Mexico through June 2006 
(lAdelman-McCarthv et al.l l2008h . WD-M star binaries are a 
byproduct of this survey which was focused on quasi- stellar 
objects. The most significant part of the analysis for this pa- 
per was completed using the SDSS spectra, covering the spec- 
tral range from 3 800 to 9 200 A at an inverse resolution of 
/1/A/l = R ^ 1800. Targets for SDSS spectroscopic follow-up 
observations were selected on the basis of photometric results. 

The spectroscopic dataset used in this study consisted of 
857 objects that were labeled as WD-M dwarf binaries during 
a search for hot subdwarf stars in the SDSS dataset. Finding hot 
subdwarfs is relatively easy, because they can be clearly distin- 
guished from normal stars by their UV-excess, a property that 
they have in common only with WDs. While the hot subdwarfs 
dominated the popula tion of UV-bright stars in earlier studies 
such as the PG survey (iGreen et al.lll986h . the depth of the SDSS 
ensures that the numbers of hot subdwarfs taper out as the vol- 
ume surveyed extends beyond the Galactic disk. White dwarf 
stars, on the other hand, are some five magnitudes fainter, and 
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their spatial proximity allows their numbers to increase steadily 
with increasing magnitude until they outnumber the hot subd- 
warfs by 100 to 1 for g > 19. Kleinman et al. (2004) presented 
a catalog of spectroscopically identified WDs in the DRl. In our 
latest detection survey of UV-bright stars in the SDSS, we de- 
tected a total of 639 hot subdwarfs and 1 1 752 white dwarf stars. 

The SDSS spectroscopic survey targets a range of objects 
based on its properties measured in the photometric ugriz survey, 
i.e. depending on whether a source is extended or point-like and 
its position in color-color space. The survey was optimized for 
detecting quasars and faint galaxies, but since these are scattered 
over significant parts of color-color space, the main strategy of 
the SDSS selection procedure is to avoid the main stellar loci. It 
is also worth noting that this procedure became more eff'ective 
after the early data release (EDR) of the survey, so the selection 
biases are not constant even within the original SDSS dataset. 

Our selection procedure resembles that of the SDSS' own 
in that we investigated a large fraction of color-color space in 
the EDR release, and refined our search area for the data sets 
of subsequent releases. Of particular importance for the WD-M 
sample was the cutoff' in g - r. For the EDR, we imposed no 
cutoff', but in later releases we imposed g - r < 0.2 in addition to 
the u - g < 0.8 criterion used to select stars with a UV excess. 
While ^ - r < 0.2 includes eff'ectively all hot subdwarf stars, 
including systems with F- to G-class companions, this limit cuts 
straight through the locus of WD-dM stars as illustrated in Fig. 
[U The truncated region of color-color space is also where the 
data for the most significant part of the quasars is found, and the 
SDSS survey therefore imposes a positive bias on WD binaries 
located in this region, while giving lower priority to objects with 
negative g - r, where single WD stars are found. However, due 
to our selection procedure, only a fraction of the WD binaries 
with g - r < 0.2 is included in this study. Also, the cool end of 
the WD population is found to be located redward of the corner 
of our selection region at u - g = 0.8 and g - r = 0.2. At this 
point, the WD sequence reaches into the population of normal 
halo stars, which the SDSS selection procedure - as well as our 
color-color cuts - avoids. A more detailed explanation of the 
c olor-color loci of WD- M star systems within the SDSS is given 
in lSmolcic etaP (|2004 . 

Our sample is made up of 857 objects, which represents 
our master sample (see Fig. [T]). We then discard noisy spec- 
tra and those in which either the WD or the red companion 
have rather weak features. Only spectra showing either a WD 
with a hydrogen-dominated atmosphere (DA white dwarfs) or 
a WD with clear helium features (DO or DB white dwarfs) 
are retained. Furthermore, we reject spectra with pollution from 
nearby light sources, which we verified with the aid of the pho- 
tometric data. WD-M pairs with cool WDs and dominant red 
components are discarded preferentially from our sample by the 
selection method, in addition to systems with significant light 
contribution due to mass overffow. The identification of our sub- 
sample of 636 systems is therefore subject to severe selection 
eff'ects. 

2.1. Premises 

As a basis for our investigations, we rely on three premises: 

1. Due to our model grid resolution of 0.5 dex for the surface 
gravity of the M star (^m) and the low resolution of the SDSS 
spectra, the uncertainties in our fits for are too large to pro- 
vide useful constraints on the secondary masses. We assume 
instead a mass-radius relation for unevolved main sequence 



stars b y using the evolutionary models of lChabrier & Baraff^ 
([T997h for a fixed M star age of 10^° yr to deduce the radius 
(Rm) and mass (Mm) of the M star from its fitted eff'ective 
temperature (T^) and metallicity [Fe/H]M. For the mass range 
considered here, evolution eff'ects on the main sequence are 
negligible within a Hubble time, thus only pre-main sequence 
stars with ages < 5 • 10^ yr might introduce larger errors in the 
radius estimate. Such young ages, however, can be excluded 
due to the cooling ages of the white dwarfs in almost all cases. 

2. Similarly, the errors in the surface gravity of the WD (^wd) 
from our fits are in most cases inadequately high for a mean- 
ingful deductio n of the WD radii. However, as shown by 
iHu et al.1 (l2007l) . most of the SDSS DAs have masses that 
cluster closely around the peak of the field white dwarf mass 
distribution at 0.58 Mq. There are observational indications 
of higher WD masses in mag netic CVs (mCVs) compared 
with t he masses of field WDs (iRamsa ^ 120001: iRitter&Kolbl 
il998i:ICr"opper et al.iil998. and references therein). The max- 
imum of the WD mass function in those binaries, although 
not as distinctive as the peak in field- WD masses, appears to 
be located between 0.7 M© and 0.8 M©. Although we can- 
not assess the orbital separations of the duets in our sample, 
except for some known close binaries (PCEBs) and optically 
resolved duets, we find no indications, such as accretion disk 
features, for a single mCV in our sample (see Sect. [6]). We 
therefore use a fixed WD mass (Mwd) of 0.6 Mq and the WD 
temperature (T^^) derived from the fit to estimate the radius 
of the WD (R wn) from the evolutionary tracks calculated by 
IWoodI (11994 . These models are available for 0.4, 0.5, 0.6, 
0.7, 0.8, and 1.0 M©. 

3. The probability of an alignment of a WD and an M star by 
chance is found to be negligible. Our reasoning is based on 
the fact that our binaries found their way into the sample due 
to the high ffux contribution of the WD to the blue part of 
the spectrum. Hence, WDs are a preferred byproduct of the 
SDSS and some have M star companions. We compute the 
probability Pu of finding at least one M star within a circle of 
diameter equ al to that of an SDSS fiber on the celestial plane, 
which is 3 " (lAdelman-McCarthv et ani2008h . by chance. It 
is given by the Poisson probability of finding y objects 

Puiv) = -e-' (1) 
v! 

and, in our case, 

Pm(v> 1) = \ -e-' , yu = A-Pm (2) 

where A = tt • (1.5 ''Y is the probed area and pu is the area 
density of M stars on the celestial plane. To estimate Pm, we 
refer to private communication with J. Bochanski, who mea- 
sured the field luminosity function of stars in the DR6. From 
there, we derived star counts in the SDSS u-, g-, r-, i- and z- 
filters for the magnitude range of our sample: 15 <u, g, r, i, z 
< 20.5. These counts yield Pm = 9.55 • 10"^ CO"^. and we get 
Pyi = 6.73 • 10"^ and an expectation value of Pm • 636 ^ 4.3 
to be the number of WD-M binaries aligned by chance within 
our sample. This is why we may use physical interaction of 
the binary constituents as a basic principle. 



3. Models 

The central element of our analysis is the synthesis of binary 
spectra from of theoretical spectra based on model atmosphere 



4 



Heller et aL: WD-M star binaries from the Sloan Digital Sky Survey 



calculations for white dwarfs and main sequence stars, respec- 
tively. Models for the full parameter range of interest were pre- 
computed, providing two grids for the spaces of white dwarf 
parameters and gwD and M dwarf parameters T^, gu and 
[Fe/H]M, respectively. White dwarf spectra were calculated for 
pure hydrogen atmospheres covering surface gravities of 7 < 
log(gwD) < 9 with a step size of 0.5 dex in a temperature range 
from 6 to 90 kK. We adopted steps of 1 kK between 6 and 30 kK, 
steps of 2 kK between 30 and 50 kK, and steps of 5 kK for tem- 
peratures above 50 kK. This core model grid was supplemented 
with existing models for extremely hot (T^^ > 90 kK) or lower- 
gravity (log(gwD) ^ {5, 6}) atmospheres and for hot helium-rich 
atmospheres. The latter were available for log(gwD) = 7.5 in 
the range of 40 kK < T^J" < 80 kK every 5kK. The grid of 
the MS spectra was almost complete for 2 600K < ev- 
ery 200 K, [Fe/H]M e {-1.5, -1, -0.5, 0, 0.3} and logCgwo) e 
{2, 3, 4, 4.5, 5, 5.5}. 

The white dwarf models were computed by D. Koester using 
his codes developed for static, plane-parallel, stellar atmospheres 
in radiative, hydrostatic, and lo cal thermodynamic eg uilibrium 
(LTE), as described in detail in iFinley et al I (1 19971) . The main 
body of this grid covers fully-blanketed pure hydrogen (DA) at- 
mospheres. These models consider convective flux according to 
the mixing length approxim ation (MLT), by using a variation of 
the standard for mulation of iMihalasI (Il978l) designated as ML2 
in the notation of lFinley et al.l(ll997r and references therein), and 
adopting a mixing length of a = 0.6 (in units of pressure scale 
height). This setup was demonstrated to be one of the most efl'ec- 
tive available config urations for reproducing DA spectra at hig h 
resolution and S/N (iHomeier et al.lll998l : iKoester et al.ll200ll) . 
This grid was extended by an existing set of hydrogen atmo- 
sphere models based on the same input physics to cover the 
more extreme parts of the white dwarf parameter space in high 
and lower gravity (sdB-like atmospheres). In addition, for 
those hot white dwarfs that could not be reproduced accurately 
by DA spectra the sequence of helium-rich (DO) models was 
used. These models assume a helium-to-hydrogen mixing ratio 
of 100:1 and the M Ll MLT version with a = 1.0 as described in 
[Jordan et"aD (1 19971) . These latter models only cover the spectrum 
up to /I < 8 000 A, and thus an important part of the dM flux had 
to be masked out in the analysis. However, none of the primaries 
were identified to be in the low-gravity or ultra-hot domain of 
the DAs/sdBs, and a more quantitative analysis of the helium- 
rich white dwarfs was beyond the scope of this work, as detailed 
below. 

The secondary spectra were calculated with version 
14.2 of the multi-pu rpose stellar atmosphere code PHOENIX 
(iHauschildt & Baronlll999) for ID spherically symmetric, static 
atmospheres in LTE, which consid er convec t ive ins tability in the 
framework of MLT according to IMihalasI (Il978l) with a mix- 
ing length parameter a = 2.0. Our mod els follow the gen- 
eral s etup used for the first GAIA grid (iBrott & Hauschildi 
120051). which included a nu mber of updates from the NextGen 
grid (Hauschildt et"aDl 19991) and the microphysics described by 
lAUard et al.l (l200lh . but ignores the eff'ects of condensate for- 
mation. We therefore restricted the range of validity of these 
models to Teff > 2800 K. Ma j or mo difications relative to the 
models of iBrott & Hauschildt (l2005l) included updating of the 
(scale d) elemental abundan ces to the revised solar composi- 
tion of lAsplund et al.l (l2005l) and considering microturbulent line 
broadening with a statistical velocity;^ = 1.0km s~\ which pro- 
vided a mor e reliable description of the observe d line profiles in 
cooler stars (iBean et al.ll2006l:ISousa et al.ll2008l) . 



4. Mathematical Treatment 

To decompose a composite spectrum, one usually fits a WD 
model to its hot region, i.e. the relatively blue wavelength range 
of the spectrum, where contamination of the companion star's 
light is rather weak. In the next step, this best-fit model for 
the WD is subtracted from the combined spectrum to obtain a 
spectrum most closely resembling the light of the secondary. 
Depending on the density of the model grid and the success 
of the preceding fit, one improves the quality of the replicated 
spectrum with each of those iteration steps, while is reduced, 
approaching the local minimum in the respective parameter 
space. Although this proceeds rapidly, there is the danger of be- 
ing trapped inside a local minimum instead of the global mini- 
mum. This procedure also implies a mutual dependence of the 
computed flux scaling factors. We therefore propose an alterna- 
tive technique. 

With our;^^ fitting method, we operate in a five-dimensional 
parameter space, spanned by T^, gwD, ^m, and [Fe/H]M. 
Each observed flux data point Ff^ in a binary spectrum with a 
total number of m observed data points is reproduced by a com- 
bination of two single- star models Xi and yt. These are weighted 
with scaling factors a and b, respectively, depending on the dis- 
tances of the stars. With the definition of from iPress et al] 
(I1992L Eq. 15.1.5), we have 



^ Z7obs z7'mod\2 / zr-obs 



X 



/ JT-obs Z7mod\^ 



(3) 



with cr, as the observational error given in the SDSS .fits file of 
the object. In the following, we use [w,] as an abbreviation for 
Wi/crj. We rewrite Eq. © as 

■■= [(Ff^-ax>-byd'] . (3) 

Since we attempt to identify that model combination of models 
X and y that corresponds to the global minimum, we have the 
boundary conditions 

= Vx\a,b) = {J^x'(a,b), ■^X\a,b)^ and (4) 

< Ax'(a,b) . 
With 



Vx'ia, b) = -2 {[xi (Ff -axi- byd] , [j,- (Ff ^ - ax, - byd]) 

(5) 



we obtain 



[F*x,]\ _ l[xj][x.y.]\la 



(6) 



which is equivalent to 



= Z" 



[Ff^xi] 
[Ff'yi] 



(7) 



This finally leads us to 



a\ _ ([Ff^Xi][y^]-[xiyi][Ff^yi] 



b) - \[Ff^y>] [xf] - [Ff^x.]j [xj] [yj] - [x. j,] 



(8) 
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Fig. 2. Example for our fi tting procedure ; Top: O bserved SDSS spectrum of J 1009 +0244 and the DA-dM combination in our fitting 
routine. As mentioned in ISilvestri et alJ (l2006h , the feature close to 5 600 A is an artifact of the SDSS data reductions, and not a 
stellar feature. Center. Decomposition of the observed spectrum. Bottom: Estimate of the fit as defined in Eq. ©. Note the mismatch 
of the emission in the line ! 



Equation ([8]) provides the scaling factors for both the WD and 
the M star model, which we then use to compute with Eq. (|3]). 
Using this procedure, we are able to avoid a mutual dependence 
of a and b since the system of equations can be solved uniquely, 
and we also avoid the eff'ects of identifying a local, and not the 
global, minimum. However, the technique is more time con- 
suming since we must first compute the distribution for the 
complete parameter ranges that we wish to consider. 
Additionally, we introduce S as an estimate of the fit, which is 
defined to be 



B(A) := 



(Fobs (/I) - ^mod(^))^ 



(9) 



and resembles the definition of x^ in Eq. ([3]). Since we do not 
evaluate the sum over all data points /, the estimate S(A) may be 
regarded as;^^g^(/l) or, in other words, as an individual estimate of 
each single data point /. It provides a means of applying a man- 
ual cutoff'-mask to noisy parts of the spectrum, implemented in 
our fitting routine. Hence, we achieve the decomposition of the 
combined spectrum, and an example of a typical binary spec- 
trum in our sample can be seen in Fig.O In Fig.O we show an 



example of the projection of the five-dimensional;^^ landscape 
onto the T^^-T^ plane. We use 



I 2 _ 2 
I J X min 



(10) 



^eff' ^WD, ^M, [Fc/HIm}, j e {up, down} 



given in IZhang et aP (Il986h to calculate the 1 cr-confidence in- 
tervals of ^ 68 % for our fits. The lower and upper grid step 
widths in the ith parameter dimension about the;^^ minimum are 
given by dij. For example, J^m ^own = 200 K = ^p. 

We define p to be the number of parameters, such that we 
have p = 5, and m as the number of data points in an SDSS 
spectrum. We have mainly m ^ 3830, but, in some cases, parts 
of a spectrum must be discarded, then m becomes smaller. We 
denote the number of statistically independent data points in the 
spectrum for an inverse resolution R, by 



dN = Rln 



(11) 
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Fig. 3. Projection of the landscape onto the T^^-T^ plane 
for J1009+0244; The best fit in this case was found at T^J^ = 
14 000 K and = 2 800K. 



where A 



min/max 



are the lower and upper wavelength limits of 
3 800 A and 9 200 A, yielding between 1 600 and 2000 resolu- 
tion elements per spectrum. We define n to be the number of 
degrees of freedom given by ^ = N - which is the number 
required for our mathematical evaluations. To ensure that the 
quality of the fits are comparable with each other, we compute 
the reduced via 
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Fig. 4. Geometrical explanation of the distance evaluation coef- 
ficient C. 



of the respective distances are comparable for all the spectra, we 
introduce the coefficient C. It weighs the displacement of a cer- 
tain binary system from the diagonal by the average distance of 
the system from Earth, which is J = (Jwd + <^m)/2 (see Fig.|4l): 



^ _ A^(JwD -^m)/2 _ -^d^ 



(JwD + dM)l2 



(15) 



Thus, a discrepancy between Jwd and du is more significant 
when the binary is closer, and we expect C to be close to zero 
for those cases in which our premises are justified and the spec- 
tral fit is good. The algebraic sign indicates the direction of the 
displacement. 



(12) 4.2. Optically Resolved Binaries 



We use the concept of the quantile that provides an estimation 
of how probable a computed is depending on the number of 
degrees of freedom n, and is defined to be 



- 1 r 

- 2'^/2r(|) Jo 



, (13) 



which provides the probability of being equal or smaller than 
a certain ;tL- 



In addition to the spectroscopic data, we also studied the pho- 
tometric data. Since we hardly detect the smeared Airy disks 
of the objects, we are unable to derive any morphological in- 
formation, but a possible separation of the spot centers allows 
us to estimate their projected mutual distance Jproj. For this pur- 
pose, we use the averaged distance d of the system from Earth, 
the displacement D of the stars on the SDSS images in units of 
pixels, and the given resolution g of the images in ''/pixel. We 
calculate the angular separation a in units of radians by means 
ofa = g Z)/3600 • ;7r/180 and by using geometry. 



4.1. Distances 

When 7?wD is inferred from th e fixed value of Mwd, the fitted 
TV^, and by using models from lWoodI (119941) . we derive the dis- 
tance JwD of the WD from Earth by obtaining the scaling factor 
a from our fitting method, which scales the model flux to the 
observed flux. Since a = (7?wd/<^wd)^, we have 



<^wD[pc] = 



(14) 



1 pc ^/a/RQ 

A similar procedure is applied to deduce the distance du of the 
M star. These two distances should be equal ideally in case the 
two stars form a physical pair. In a plot with Jwd on the abscissa 
and du on the ordinate, physical systems should lie on or near 
the diagonal where Jwd = du- Their displacement from the di- 
agonal is given by ^/2 (Jwd - <^m)/2. To ensure that the estimates 



Jproj = 2(itan^|j 



(16) 



The projected distance, however, is shorter than the true spatial 
separation r (see Fig. [5]). If we knew r of a system and were able 
to measure Jp^oj many times during an orbit, or if we knew r for 
many systems and were able to measure the respective Jproj for 
each system, then we would state that, on average, Jpj-oj reaches a 
value closer to 0.637 r, the more observations that we have. We 
can deduce this mathematically by considering the mean value 
of (ipj-oj, which is 



< ^proj > = < r sin(^) > = r < sin(^) > 



(17) 



where the mean of the sine of the projection angular (p is given 
by 
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line of sight 




Fig. 5. Jproj is the observed projection of the displacement r of the 
two stars. Averaged over many observations, the latter is about 
1.57 times smaller than the real spatial separation. 



< sin(^) > 



dip sin(v?) 



- ^ 0.637 

71 



(18) 



This implies, that the true spatial separation is, on average, about 
1/0.637 ^ 1.571 times longer than the observed projected dis- 
tance. We will denote this statistical correction of Jproj with 

jcorr •- 1 571 J . 
"proj • I ^ "proj. 

Using the optical separation of the stars, we can use the vis- 
ible separation and estimate the period of the system by means 
of 



5. Results 

The temperature function for our WD sample can be seen in Fig. 
[71 As mentioned above, it is subject to considerable selection 
effects. While the strong peak at 17 000 K belongs to the DAs, 
the second and smaller peak at 40 000 K, which is the lowest 
temperature in our DO model grid, stems from the cool helium- 
rich WDs. We convolved the original data distribution with a 
Gaussian probability function about each data point, to account 
for the uncertainty in the derived parameters. As a measurement 
of standard deviation, we used 2 000 K and normalized the distri- 
bution to 1 . Since most of the WDs in our sample have a surface 
gravity of between 7 and 8 dex (see Fig. [11]), the shape of the 
WD-distance function in Fig. [9] is similar to that of the tempera- 
ture function. 

In Fig. [SI we show the temperature function for the M stars. 
Although the standard deviation is merely a tenth of that for T™, 
the function is not so smooth because of the relatively low num- 
ber of model grid points. The tail towards temperatures higher 
than those of M stars is due to 4 objects at = 4 000 K with 
insecure fits, i.e. because of He emission in the spectrum or large 
distances to Earth and therefore weak spectral features. Our re- 
sults for the other free parameters, and [Fe/H]M, can be seen 
in Figs. [12] and [131 while the distance distribution for the M stars 
is given in Fig. [T0|_ 



As shown in Hu et al. I (l2007h . the standard deviation from 
the maximum of the WD mass function at 0.58 Mq is 0. 1 Mq. 
This corresponds to an error in our WD distances of 17%, and 
this uncertainty creates dispersion in the data from the diagonal 
in the dwD-du plot, namely 



P > 2n 



proj 



' G (MwD + Mm) 



(19) 



where we assume that Mwd = 0.6 Mq on the basis of argu- 
ments given above and for the sake of consistency with the fit- 
ting pro cedure. We derive M m from evolutionary tables calcu- 
lated bv lBaraff'eetaD (Il998h (see Sect.[rT]). Equation ^ may 
only provide an approximate approach, due not only to the un- 
certainties in the masses but also because of the approximation 
of the projection of the mutual distance is smaller and almost 
equal to the sum of the major axes ^wd and Gm of the ellipses on 
which the stars move: dp,oj ^ ^wd + ^^m- If we take into account 
the averaged correction of the projection, then the corrected pe- 
riod P""'^'^ is longer than P. We note that deviations dMwD from 
Mwd = 0.6 Mq have a significant eff'ect, since the mass was al- 
ready used for the derivation of Jproj. Since d^^^- oc r^^^ oc M^^, 
we derive 



dP 
dMwD 
which leads to 



dM^^ 



Mwd + M^ 



1/2 



2Mwd + Mm 

dP oc T-rdMwD 

2 [MwD(MwD + MM)f' 
-M^^ dMwD 



(20) 



(21) 



Mwd > Mv 



Thus, deviations in the period are more significant when the 
mass of the white dwarf is lower and, for Mwd = 0.6 Mq, they 
could be considerable. However, since dMwD is usually small, 
impacts due to variations in the WD mass should not have a sub- 
stantial impact on the assessment on the lower limit to P, i.e. 
approximately 10 % for P. 



C = V2 



0.17 JwD 



^WD + 

and for physical systems with Jwd ~ du 



C 



V2 — 



0.12 



(22) 



(23) 



Taking into account all the other error sources mentioned before, 
we would expect the systems to be mainly distributed within the 
fan of C = ± 0.25 to both sides of the diagonal. In Fig. [6l we 
present the derived distances for all systems. Indeed, this plot, 
and in particular the zoom-in section in the right panel, shows a 
scattered distribution, fanned out to both sides of the diagonal, 
with most of the binaries located in the C = 0.25 -tolerance inter- 
val. However, we also detect a trend towards higher distances for 
the WD component - the scattering is asymmetric with respect 
to the bisecting line. 

For all of the analyzed spectra, = I (see Eq. ([T3])), 
which implies that a equal to or smaller than the calculated 
Xmin is probable for the respective number of degrees of free- 
dom. In other words, the quality of our fits is poor in mathe- 
matical context. The standard deviations of the measured pa- 
rameters are quite weak in terms of physical significance. As 
an example, we consider J1009-h0244 (see Fig. [3]). We com- 
pute (TrWD,down ~ 80K ^T^eT'UP' <^r^.,down ~ 5K ^ ^T^eff,"?' 

which demonstrates that the systematic errors are larger than the 
mathematical ones. For a conservative estimate of the errors due 
to the incomplete molecular data in the M star models, SDSS 
flux calibration errors a nd interstellar reddening, we refer to 
iHiigelmeyer et al.l (l2006l) . and assume an uncertainty of 2 000 K 
for the WDs with temperatures smaller than 50 000 K and a l-cr 
interval of 100 K for the MS stars. For T^J" > 50 000 K, the ab- 
solute value of our accuracy is given by half of the model step 
width of 5 000 K. Our accuracy in the surface gravities is limited 



8 



Heller et aL: WD-M star binaries from the Sloan Digital Sky Survey 



by the low resolution of the SDSS spectra and the step size of our 
model grid, and we therefore have cri^gig^^) ~ 0.5 dex ^ criog(gM)El. 
The metallicity determination is only accurate to about 0.3 dex. 

5.1. White Dwarfs Showing He Lines 

From our original input sample, 616 spectra can be fitted well 
using a DA primary. However, among the remaining systems. 
He features are evident in a small number in the primary spec- 
trum. These are most readily apparent at relatively high 
since spectra with distinctive DB features were efficiently re- 
moved in the initial selection. We found that 20 of the remain- 
ing spectra, which were dominated mainly by the He ii Bracket- 
equivalent series, were fitted more accurately by spectra from 
our DO model library. These stars cluster predominantly at the 
cool end of the DO sequence: 16 of 20 are fitted with the lowest 
effective temperature in the DO grid of 40 000 K. Many of the 
them show stronger Hei lines in addition to Hen, and would 
probably be classifie d more accurate l y as D BO, according to 
the nomenclature of 'Wese mael et al.l (Il993l) . We suspect that 
this distribution does not reffect the true luminosity function of 
Helium-rich white dwarfs, but is rather biased by a contribution 
of WDs below the T^^ limit of our DO models. Since the number 
of these systems, and an additional number of potential DB sys- 
tems that could neither be fitted by DA nor DO models and were 
omitted from our catalog, is too small for a meaningful statistical 
analysis, we have not attempted to obtain more reliable fits with 
additional DB models. 

Three DOs in our sample (J075 6-F4216, J09 16-^0521, an d 
J1336-0131) were already studied bv lHiigelmever et al.l (l2006h . 
For the first one, the parameter values derived in our study are 
in good agreement with those of the cited publication. While 
they found the white dwarf to show an effective temperature of 
52.5 ± 0.8 kK at gwD = 7.4 ± 0.05 dex, we find = 50 ± 5 kK 
and gwD = 7.5 ± 0.5 dex. For the MS companion, they derived 
= 3.2 ± 0.2 kK, gM = 5 ± 0.5 dex, and Mm = 0.18 M©, 
whereas we measure = 3.4 ± 0.1 kK, = 5 ± 0.5 dex, and 
Mm = 0.27 Mo. 

The second binary , J0916-h0521, was proposed by 
iHiigelmeyer et al.l (l2006h to be one of the few optical pairs 
aligned by chance. With C = 0.461 for that system, we tend 
to confirm their findings although the data points for many 
systems are located in the same region of the Jwd-<^m diagram 
around J0916-h0521. The other parameter values do not differ 
significantly from each other, except for the M star masses, for 
which they derived 0.51 M©, while we find Mm = 0.18 M© due 
to ^ 10% lower T^^ and subsolar metallicity found in our fit. Our 
results f or the third systeni , J133 6-0131, agree with the values 
given by IHiigelmeyer et all (l2006h . 



5.2. Spectra with Emission 

In the spectra of 167 objects in our sample we find at least H^ 
emission, and some of them also exhibit other types of B aimer 
emission. Most of these systems probably harbor an active M 
dwarf with chromospheric emission, or the respective M star 
might be irradiated by the primary, and experience photoioniza- 
tion and recombination. These systems are non-accreting since 
they do not show He i or He ii features. A more detailed treatment 
of em ission in WD-M star binaries is given in Silvest ri et 
(I2006h . For an assessment of the magnetic activity, we scanned 



^ iReinersI (l2005h provided an even more conservative error estimate 
of 1 dex, due to the low resolution of the spectra. 



our fits for pronounced peaks in 5, as defined by Eq. (|9l). An ex- 
ample is shown in the lower panel of Fig.[2l Our M dwarf models 
did not consider H^ emission, and therefore the spike in the fit- 
ting error at the H« wavelength was - if present at all - always 
related to an emission feature in that line. We decided to attribute 
the peak in S close to 6 563 A to H« emission, if it exceeded the 

ambient noise level by at least a fraction of t hree. 

In their study of 101 WDMSs, .Rebassa-Mansergas et al.l 

(l2007h identified 18 PCEBs and PCEB candidates of which 12 
are included in our master sample (see Table [3]). Eleven of them 
are in our spectroscopic sample and 8 of them are marked as 
H emitters therein. Due to the strong contamination of the 
PCEB (candidate) spectra our fits are often inconsistent in terms 
of dwD and du, whereas they should ideally be equal. If the pa- 
rameter values of our best-fit solution differed significantly from 
those provided by Rebassa-Mansergas et al. ( 2007), we chose to 
adopt instead the parameters corresponding to their solutions. 

As an example, we refer to J0052-0053, which shows an 
ambiguous x^^d distribution. This is probably due to the weak 
WD and M star features in the spectrum and due to the large 
distance of the system. The results of our fitting routine for 
J0052-0053 are T^J" = 29 000 K, log gwD = 8.5, and a distance 
to Earth of 1 612 pc for the WD. The corresponding parameters 
for its companion are computed to be = 3 200 K, log ^m = 
5.0, [Fe/H]M = 0.0, and dwD = 223 pc. The respective distance 
evaluation coefficient of C = 1.07 emphasizes the grave disso- 
nance in the objects' distances from Earth. It should be close to 
zero for consistent fits as defined in Eq. ([T5l) . Hence, we fix the 
parameters for the WD to be Mwp = 1.0 Mn, = 15 000 K 
and log gwD = 8.5 - following iRebassa-Mansergas et al.l (l2007l) 
- and set = 3 600 K to enforce du ~ ^^wd- Using this pro- 
cedure, we deduce Jwd = 579.9 pc and du = 577.2 pc, in good 
agreement with the value of 505 pc for the distance of the WD 
derived by lRebassa-Mansergas et al.l (l2007l) (see Table [2] for the 
whole result). 

We also find 4 systems showing Hei and Hen emission, 
in addition to H^ features, that harbor a hot primary (T^^ > 
65 000 K). Since Hei is not typically observed in active M 
dwarfs, these objects might either be CVs, i.e. binaries that un- 
dergo mass transfer, or their mutual separation might just be 
large enough as to avoid mass overflow on the one hand but 
close enough to produce the emission lines by irradiation effects 
on the other hand. J124 9-h0 357 was actu a lly cla ssified as CV 
bv lDownes et al.l (l200lh and lSzkodv et al.l (l2004l) . However, its 
CV status can be doubted since the spectrum shows no features 
of an accretion disk and the He emission lines are not dynami- 
cally broadened by a potentially rotating disk. The emission fea- 
tures in J1317-h6731 and J1439-0106 are most likely due to ir- 
radiation of the M component. Their spectra show strong H^ as 
well as He i and He ii emission . Both objects were already men- 
tioned in lSilvestri et al.l (l2006h and the detec tion of the latter one 
was originally published in [Raymond et al.l (|^03). The system 
J2 125 -01 07 that we fitted with a DO primary from our model 
grid is actually known to be a close systern with a PG1159 pri- 
mary (iNagel et al.ll2006l: ISchuh et al.l[2008l [2009h . We misclas- 
sified that object since we did not include PG 1159 models in 
our repertoire. Our value for of 4 000 K should be taken with 
reservation due to the known strong irradiation processes from 
the hot primary on the MS companion. 

From the 167 spectra with H emission features, we subtract 
8 known PCEB candidates and the 4 systems from the previous 
paragraph, which yields 155 WD-M star binaries, corresponding 
to a fraction of 24.4 %, that probably harbor an activ e M dwarf 
Although this value matches the number found by IWest et aO 
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Fig. 7. Temperature function for the WDs; While the major peak 
at 17 000 K belongs to the DAs, the bump at 40 000 K is due to 
the preferential selection of cool DQs for our sample. T his plot 
looks very similar to that shown in lSilvestri et alJ (I2006h . except 
for the DO feature and that it is smoothed. 



(I2OO4I) for active field M dwarfs in the SDSS, our sample of 
MS stars mainly consists of stars with spectral types earlier than 
M5 (see Fig. [15]). For this range, and in p articular for dM st ars 
earlier than M4, the activ ity fraction from IWest et al.l (|2004|) is 
much smaller than 10%. ISilvestri et al.l (l2006l) detected a sig- 
nificantly higher fraction of active M dwarfs in WD-M binaries 
(between 20 % and 50 % from MO to M5), which they attributed 
to a rotational spin-up of the secondary due to gravitational in- 
teraction with the WD. They also reported a lack of active dM 
stars with ages > 0.8 Gyr. Since they imposed a less stringent 
(g - r) cutoff' (< 0.7) for most of their sample than we had (< 
0.2), hotter, and thus younger, WDs passed their selection pro- 
cess, which would increase their fraction of active M dwarfs to a 
higher value than ours. Nevertheless, our results for the fraction 
of active M dwar f s is in better agreement with the results found 
bylSilvestri et al. (l2006h than it fits to the low fraction found by 
IWest etal.1 (12004 . Even if we consider only unresolved binaries, 
as lSilvestri et al.l (l2006h did, our fraction of binary systems with 
emission that probably originates in the M component does 
not decrease below 23.2 %. 



5.3. Optically Resolved Binaries 

Of our 636 binaries, 41 were chosen for follow-up studies be- 
cause, firstly, the red star was found to be located within an area 
around the WD that was covered by the SDSS fiber and, sec- 
ondly, these stars had a sufficiently wide separation to enable us 
to distinguish between the two components. Figure fT4l provides 
a typical illustration of these systems. For photometric binaries 
with a separation clearly larger than the SDSS fiber radius of 
1.5 '\ we assume that the M star on the SDSS image is not the 
one represented in the spectrum. The respective fitting results 
are indicated in Table |2] and should be interpreted with caution 
since we may be unable to determine if a significant fraction of 
the light contribution from the red star was collected by the fiber, 
despite the separation being larger than 1.5 The results for the 
41 clearly separated objects can be seen in Table [T] The typical 
projected distances are large and of the order of some hundred 
AU with a mean value of roughly 650 AU. The widest separa- 




Fig. 8. function; The cutoff' at 2 600 K towards lower temper- 
atures is due to our lack of cooler models. The peak at 3 200 K 
is probably generated by SDSS selection eff'ects and does not re- 
ffect the true M dwarf population. The tail towards > 3 800 K 
is caused by 4 objects with insecure fft at = 4 000 K and the 
Gaussian decay. 
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Fig. 9. The distribution of the WD distances in our sample is a 
consequence of the restricted magnitude range of the SDSS. It 
has its maximum at 354 pc. 



tion that we find is 1 700 AU for J1006-F5633, where the true 
spatial separation should be even wider. The orbital periods that 
we derive are typically of the order of 10^ yr with J1006-h5633 
indicating a period of > 70 000 yr. 

5.4. Very Low-Mass Objects 

We ffagged the secondaries in our sample to be very low-mass 
objects (VLMOs), provided that they had a mass < 0.1 Mq and 
logigu) = 5, or Mm < 0.09 M© and logigu) > 4.5, to iden- 
tify potential candidates for substellar companions. Seventy M 
stars in our reservoir are VLMO according to that definition (see 
the electronic version of this paper for the table of results). The 
M masses were calculated from the fitted and [Fc/H Im us- 
ing evolutionary tables calculated by lBaraff'e et al.l (1 19981) . Their 
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Fig. 6. Distribution of the derived distances of the stellar components of our systems; The diagonal is the ideal for physical binaries. 
The dotted curves span a tolerance fan for C = 0.25. Left: Note that for very far away - and thus rather faint - objects with a distance 
of more than ^ 1 000 pc the distance for the respective binary mate is typically very different. This is due to the weak spectroscopic 
features of the distant stars and the low spectral resolution. Right: Same as left but zoomed in. 
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Fig. 10. The distribution of the M star distances in our sample 
is a consequence of the restricted magnitude range of the SDSS. 
Its maximum is at 334 pc. 
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Fig. 11. Distribution of the WD surface gravities. 

models would be expected for these objects, biasing the radius 
estimate and thus also the assessment of the distance. 



models are available for [M/H] = -2, -1.5,-1.3,-1, -0.5, and 
0, while no specific tracks exist for those M stars, which are 
fitted to have [M/H] =0.3. Since radius and mass at constant 

decrease monotonously with metallicity in the other tracks, 
we can only quote the [Fe/H]M = values as a lower limit for 
Mm and Rm, and consequently also for du- However, there is 
only one VLMO candidate in our sample that has a metallicity 
of +0.3: J1323+3018. The other 69 VLMO candidates are most 
accurately fitted with [M/H] < 0. 

If any of the VLMOs were found to be substellar, they would 
of course by definition be in the cooling phase of their evolu- 
tion, and more significant deviations from the old main-sequence 



6. Discussion 

The M star temperatures that we fitt ed are in good agreement 
with the spectral types derived by iRebassa-Mansergas et al.l 
(2007). In Fig. [151 we compare our fits for with their spec- 
tral types for 41 objects (left panel ) and with the spectral classi- 
fication from lSilvestri et al.l (l2006h (right panel) for 446 objects, 
which were both in our own and their sample. The authors of the 
latter publication derived their dM types on t he basis of t e mplat e 
spectra and color indices as described in H awley et al.l (l2002h . 
Since we project the number of M dwarfs per and per spec- 
tral type onto the plane spanned by these two parameters, this 
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Fig. 12. Distribution of the M star surface gravities; The few ob- 
jects with < 4 dex are probably M giants. 



CT= 0.3 




Fig. 13. Metallicities of the M stars in our sample, normalized 
to 1; The plot shows a distinct peak at [Fe/H]M = 0. PHOENIX 
models for [Fe/H]M > 0.3 were not available. Among the 70 
candidates for a VLMO, there is only Jl 323+301 8 showing a 
metallicity > 0. 



plot is related to Fig. [8l In both of the figures, the maximum is 
reached at (3 200 K , M4), and the plot shows a monotone de- 
crease of the spectral type with increasing temperature. These 
counts probably do not relate directly to the true or spec- 
tral type function due to selection eff'ects such as the increasing 
number of M dwarfs towards later types on the one hand, and 
decreasing visibility of the secondary component on the other 
hand. 

In the iReb assa-Mansergas et al.l (l2007|) study, the MS stars 
showed a tendency to be at larger distances from Earth than the 
WD primary (see Fig. O. The authors ascribed this trend to 
spots on the secondaries' surfaces, which would cause the spec- 
tral types to appear too early. The radii derived from the spectral 
type-radius relation would then have been too large, in addition 
to the computed distances of the MS stars. In our study, we ob- 
serve the opposite trend, namely du being smaller than Jwd, and 
we consider three possible eff'ects contributing to that trend: 



Fig. 14. As an example for one of the 41 optically resolved 
binaries with further treatment, we show the SDSS image of 
Jl 127-0028. The image resolution is 0.025 'Vpixel. In this 
grayscales image, the red companion is set black while the larger 
and brighter spot is the white dwarf. 



1 . iReinersI (l2005l) showed that the absorption in the TiO 6-band 
is systematically underestimated by PHOENIX M star spectra. 
Since as a primary temperature indicator, this band system 
deepens with decreasing T^, model spectra fits to this fea- 
ture may also underestimate systematically T^. This eff'ect, 
if present in our models as well, would cause du to be sys- 
tematically underestimated. Figure [2] shows indeed that the 
model predicts a relatively strong TiO ^-band at 8 870 A com- 
pared with the 6-band at 8 450 A. However, the mismatch is 
not severe, and the y-band at 7 055 A is reproduced well. We 
therefore do not expect uncertainties in the molecular opaci- 
ties to introduce a strong bias in our results. 

2. As mentioned in Sect. 12.11 there is observational evidence 
that the masses of WDs in magnetic, accreting binaries 
are mainly between 0.7 Mq and 0.8 Mq, instead of around 

0. 6 Mq as for field WDs. There is no evidence for accre- 
tion in any of the spectra in our sample. Only 1 1 out of the 
636 objects presented here are known to be close systems, 

1. e. PCEBs, while 41 are quite widely separated as inferred 
from the SDSS images. We cannot assess of the remaining 
orbits on the basis of our data, and we cannot state if the WD 
masses are close to 0.7 or even 0.8 M©. However, a trend to- 
wards higher WD masses would be compatible with a bias 
towards 

^wD ^ dy[. An increase in mass leads to a decrease in 
radius for a WD, and consequently a smaller distance to Earth 
is required to reproduce the observed flux. Our possible un- 
derestimation of MwD could thus contribute to a systematic 
overestimation of Jwd- 

3. Based on the efl'ective temperature and the metallicity from 
our spectral fits, th e M sta r radii that we deduced from the 
IChabrier & Baraff'd (1 19971) model tracks are systematically 
underestimated for low-mass objects (Mm < 0.3 M©, see 
Fig. [16]). This discrepancy between theory and observations 
is well known in the research field of low-mass stars and 
emerged in observati ons of eclipsing b inary systems with a 
low-mass component (iRibas et al .1120071) . With our - in statis- 
tical terms - large sample, we support this claim. We encour- 
age the reader to look at Fig.[T6]to see that the trend vanishes 
for MS stars with masses 0.3 MqI 

The VLMOs in our binary systems (Sect. 15.41) are auspicious 
targets for follow-up investigations of their mass with time- 
resolved spectroscopy. Unambiguous mass determinations, in- 
dependent from the orbital inclination of the system, are possi- 
ble in eclipsing binaries, where the inclination can be taken from 
the eclipse light curve or can be assumed to be close to 90 °. For 
systems whose orbits can be measured by astrometry, the incli- 
nation might also be more tightly constrained; although 4 of the 
VLMOs presented here are also optically resolved (see Table 
[B, time-resolved astrometric measurements would be unhelpful 
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Fig. 15. Left: Comparison of the spectral types given in lRebassa-Mansergas et alJ (l2007|) with the effective temperatures we deduced 
for 41 M stars, that match both their sample and ours; 11 of them are marked as PCEBs, as confirmed i n their study via RV 
measurements. The labels indicate the counts per grid point. Right: Comparison of the spectral types given in lSilvestri et al.l (l2006|) 
with the eff'ective temperatures we derived for 446 M stars, that match both their sample and ours; The number of M dwarf stars per 
grid point is projected onto the plane spanned by and the spectral type. The outer contour marks the path of three counts and 
each subsequent contour symbolizes an increase of three counts. The maximum is at (3 200 K , M4) with 115 counts. Their average 
spectral mismatch was ±1, while our 1-cr accuracy for is 100 K. 
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Fig. 16. Left: Distribution of the derived distances of the stellar components from lRebassa-Mansergas et al.l (l2007l) : In contrast to 
our results, their binaries tend to have a larger value for du than for Jwd- Right: For M stars with masses > 0.3 M© there is no 
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since the orbital periods are of the order of 5 000 yr. Furthermore, 
even using high-resolution astrometry with adaptive optics, it 
would be difficult to resolve systems of less than 10 yr orbital pe- 
riods, since virtually all of our targets are at distances > 100 pc. 

There is a relatively large number of widely- separated bina- 
ries with orbital distances > 250 AU and up to 1 700 AU: 41/636 
= 6.4%. This is, however, not a stable assessment; we suggest 



five eff'ects that may smear the optical-binary fraction to either 
higher or lower values. 

1. If two objects are so close together that they appear as one 
elongated object, the SDSS classification procedure may pre- 
select it as a possible galaxy. This would imply that we have 
underestimated the true optical-binary fraction. 

2. A small number of physical pairs aligned along the line of 
sight could also drive the true optical-binary fraction to a 
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higher value. We would have failed them to detect because 
of their unfavorable geometrical constellation. 

3. Some of the stellar duets on the SDSS images with angular 
distances > 1.5'' could also be physical pairs (but the sec- 
ondary on the image would not be the one represented in the 
respective SDSS spectrum). 

4. A significant contingent of binaries with mass overflow from 
the MS companion to the WD are supposed to be located out 
of the color-color region that we studied in compiling our 
sample. These stellar duets are close systems with no optical 
separation on the SDSS images, which pushes the optical- 
binary fraction to lower values. 

5. The statistical considerations presented in Sect. [2?T] demon- 
strated that about 4 of the binaries in our sample should be 
aligned by chance - without a common evolutionary back- 
ground. Probably none, or at most one, of the optical binaries 
is one of these outliers. 

To help consider whether stellar duets with such large orbital 
separations are stable over long times cales, we refer the reader 
to the paper of [Weinberg et al.l (1 19871) and Fig. 2 therein. Their 
calculations, including both stars passing by and encounters with 
subclumps within giant molecular clouds, showed that bina- 
ries with a total mass of 1 Mq and initial orbital separations of 
around 650 AU have a typical lifetime of more than the age of 
the Universe. 
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Table 1. Optically resolved binaries (41 objects); We use the mathematically exact values of Jproj and Jp°^j to compute the lower limit 
for P and the statistically corrected F'''^\ But for this table, distances greater than 500 pc or 500 AU are rounded to the next hundred 
and else to the next fifty. Periods are rounded to the next 500 yr for P < 10 000 yr, to the next 1 000 yr for 10 000 yr < P < 50 000 yr 
and else to the next 5 000 yr. ^: VLMO candidates. 



Designation 


^proj [ ] 


d [pc] 


Jproj[AU] d^Z'^AV] 


Mm [Mo] 




pcorr [yj.] 


J0017-0009 


1.4 


500 


700 


1000 


0.46 


16000 


33 000 


J0122+1542 


1.0 


300 


300 


500 


0.17 


6500 


13 000 


J015 1-0800 


1.0 


350 


350 


600 


0.17 


8 000 


16000 


J0215+1418 


0.9 


700 


600 


900 


0.27 


16000 


30000 


J0249+3342 


1.4 


500 


700 


1200 


0.17 


23 000 


45 000 


J0348-0614 


1.1 


500 


600 


900 


0.17 


16000 


31000 


J0725+4145 


1.3 


400 


500 


900 


0.17 


15 000 


29 000 


J0729+4304 


1.5 


150 


250 


400 


0.12 


5 000 


10000 


J0739+2743 


0.6 


700 


350 


600 


0.17 


8 000 


16000 


J0740+3859 


1.1 


300 


350 


600 


0.12 


8 000 


15 000 


J0741+3808 


0.9 


300 


250 


400 


0.17 


5 000 


9 000 


J0752+4332 


1.3 


600 


800 


1200 


0.18 


23 000 


46000 


J0800+5002 


1.1 


700 


800 


1200 


0.27 


22 000 


44000 


J0801+2216 


1.3 


500 


700 


1 100 


0.17 


20000 


40000 


J0806+4035 


1.2 


350 


450 


700 


0.17 


11000 


21000 


J0809+1251 


1.3 


400 


500 


800 


0.17 


14000 


27 000 


J0813+2152 


0.9 


900 


800 


1300 


0.18 


28 000 


55 000 


J0829+2701 


1.4 


1000 


1400 


2200 


0.27 


55 000 


110000 


J0845+2348 


1.1 


1400 


1400 


2 200 


0.46 


50000 


100000 


10904+5621^ 


1.5 


200 


300 


450 


0.10 


6000 


12000 


J0931+3941 


1.4 


350 


500 


800 


0.17 


13 000 


26000 


J0939+5729 


1.3 


250 


300 


500 


0.17 


6000 


12 000 


J0942+1846 


1.3 


600 


700 


1200 


0.17 


23 000 


45 000 


JlOOl+3203 


1.3 


250 


350 


500 


0.17 


7 000 


13 000 


J1006+5633 


1.4 


1200 


1700 


2700 


0.46 


70000 


135 000 


J1032+3722 


1.3 


350 


450 


700 


0.17 


10000 


20000 


11127-0028''^ 


1.5 


200 


250 


400 


0.10 


4 500 


9 000 


Jl 127+4249 


1.2 


350 


400 


600 


0.12 


9 000 


18 000 


J1205+0312 


1.4 


250 


400 


600 


0.17 


8 500 


17 000 


11209+6510^ 


1.4 


250 


350 


500 


0.10 


7 500 


14 000 


J1210+0549 


1.4 


350 


500 


800 


0.17 


12 000 


24 000 


J1216+4328 


1.5 


900 


1300 


2 000 


0.46 


46000 


90000 


J1242+4506 


1.5 


450 


700 


1 100 


0.17 


20000 


40000 


J1253+5813 


1.4 


400 


600 


900 


0.27 


14000 


28 000 


J1304+1449 


1.3 


500 


600 


1000 


0.17 


18 000 


35 000 


J1347+4129 


1.3 


1200 


1600 


2 500 


0.46 


65 000 


125 000 


J1456+4824 


1.4 


450 


600 


1000 


0.17 


18 000 


35 000 


J1606+4217 


1.1 


500 


600 


900 


0.17 


17 000 


35 000 


J1630+1302 


1.4 


600 


900 


1420 


0.17 


31000 


60000 


J1744+2442 


1.4 


1 100 


1500 


2400 


0.46 


55 000 


115 000 


12200-0715^ 


1.5 


200 


250 


400 


0.10 


5 000 


10000 



le 2. Example table of results for 636 WD-M star binaries; Typical uncertainties are mentioned in the text. The 1 1 PCEBs and PCEB candidates from lRebassa-Mansergas et ID 
i?!), as well as the two CV candidates mentioned in Sect. 15.21 are shaded in gr ay. * : M component probably a g iant due to the very low surface gravity of 3 dex. ^: VLMO 
iidates. H\ and other Balmer line emission, rm*: WD parameters taken from lRebassa-Mansergas et aP (l2007l) . •: Clearly resolved objects with projected angular distance 
5 " . o: Objects with slight angular separation and without further treatment. *: Objects separated by more than 1.5 '\ likely to be not those represented in the spectrum. ?: The 
ometric data corresponding to the object id specified does not exist in the SDSS database. The complete table is available as Supplementary Material to the online version of 
article at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsweb.u-strasbg.fr/cgi-bin/qcat7J/A-hA 



^ignation 


JwD [pc] 


dM [pc] 


C 


[K] 


Cff [K] 


^wd[^o] 


Rm[Rq\ 


log gWD 


log gM 


Mwd[Mo] 


Mm[Mo] 


[Fe/H]M 


Emission 


WD 


Res. 


m 


N 


red,min 


17-0009 


475.8 


514.6 


-0.055 


50000.0 


3600.0 


0.016 


0.433 


7.5 


4.5 


0.6 


0.46 


-0.0 


H 


DA 


• 


3783 


1572 


4.29 


17+0040 


547.2 


253.5 


0.518 


13000.0 


3200.0 


0.013 


0.189 


7.5 


5.0 


0.6 


0.17 


-0.0 




DA 


* 


3824 


1599 


4.10 


39+2548 


370.9 


335.9 


0.070 


34000.0 


3200.0 


0.015 


0.189 


8.0 


5.0 


0.6 


0.17 


-0.0 




DA 




3838 


1599 


2.26 


41+1511A* 


1332.0 


607.1 


0.529 


32000.0 


3600.0 


0.014 


0.190 


7.5 


3.0 


0.6 


0.18 


-1.0 




DA 




3842 


1602 


2.74 


52-0053* 


579.9 


> 577.2 


< 0.003 


15000.0 


3600.0 


0.008 


> 0.433 


8.5 


3.0 


1.0 


>0.46 


+0.3 


//,rm* 


DA 




3843 


1596 


19.69 


54-0025 


385.1 


> 383.1 


< 0.003 


17000.0 


3200.0 


0.015 


> 0.189 


8.0 


4.0 


0.5 


>0.17 


+0.3 


H, rm* 


DA 




3843 


1596 


6.51 



ie 3. Excerpt of our catal og, the master sam ple, of 857 WD-M binary stars , s: ISilvestri et al.l (l2006l). r : iRavmond et all (|2003|). kl : iKleinman et al.] (12004)), m : 
^ssa-Mansergas et al.l (l2007h (non-PCEBs), rm*: iRebassa-Mansergas et aD (l2007h (PCEB), e: lEisenstein etaP (l2006h . h: iHiigelmever et al.l (l2006h . g: iGreensteinI (Il986h . 
here was no alternative name found. The complete table is available as Supplementary Material to the online version of this article at the CDS via anonymous ftp to 
rc.u-strasbg.fr (130.79.128.5) or via http://cdsweb.u-strasbg.fr/cgi-bin/qcat 7J/A+A 



)rt Name (SDSS) Full Name (SDSS) 



Spectrum File 



References alt. Names (as found in the SIMBAD database) 



•01+0006 
•17+0040 
•17-0009 
•26+1444 
•36+0700 
•39+2548 



J000152. 10+000644.7 
J001733.59+004030.4 
J001749.25-000955.4 
J002620.41 + 144409.5 
J003602.59+070047.3 
J003925.22+254823.7 



spSpec-51791-0387-157.fit 
spSpec-51795-0389-614.fit 
spSpec-51795-0389-112.fit 
spSpec-52233-0753-079.fit 
spSpec-53709-2312-164.fit 
spSpec-53327-2038-380.fit 



s6, s7,r,kl SDSS J001733.59+004030.4 (Oa.'O) 

s6, s7,r,kl,rm SDSS J001749.25-000955.4 (Oa' 1) PB 5848 (Hr5) 

s6, s7 PHL 2888 

PB 6052 PB 6052 (08?!) 

X 



